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Abstract 

We investigate the decay and the production mechanism of the resonance X(1812) recently 
observed in the J/tp -> (1812), X(1812) -> u<f> at BESII. The decay widths of X(1812) -> 
7777', r/77, u>4>, K + K~, p + p~, wll>, K* + K*~ and tt + it~ are evaluated based on the scenario of the 
X(1812) as a candidate of (K*K*) molecule. It turns out that the quark exchange mechanism plays 
an important role in the understanding of the large decay width for the X(1812) — + uxfi. It is also 
found that the decay widths for X(1812) — > 7777' and r\r\ are enhanced by the quark exchange mech- 
anism. These channels are suggested to be the tools to test the molecular scenario in experiment. 
The branching fraction of Br(X — > u>4>) is evaluated to be about 4.60%. Searches for additional 
evidence about the X(1812) in J/tp radiative decays are reviewed. In the molecular scenario, the 
X(1812) production rate is also evaluated to be T(J/t/j -> jX)/T(J/ip -> jK*+K*-) = 2.13l^, 
which is close to the measured value 2.83 ± 0.92. 

PACS:13.20.Gd, 12.39.-x, 13.30.Eg. 
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1 Introduction 



The observation of a near-threshold enhancement in J /%p — > jX, X — > vcj) at BESII pQ immediately 
provokes discussion about its nature. This enhancement is reported to favor J PC = ++ in a partial 
wave analysis with a mass and width of M = 1812±^(stat) ± 18(syst) and T = 105 ± 20(stat) ± 
28(syst) McV, respectively, and a production branching ratio, B(J/ip — * jX) ■ B(X — * cj(f>) — 
[2.61 ± 0.27(stat) ± 0.65(syst)] x 10 -4 . This resonant state is named as JT(1812) in the following 
part. 

Emergence of the X(1812) adds a new puzzle to the scalar sector of mesons. The most interesting 
feature about X(1812) is that it seems to have a strange production and decay mechanism observed 
in the J/ip — > jujip. If the valence quarks of ui and cf> are assigned as (uu + dd) / y/2 and ss, one 
expects that this double OZI suppressed process should have a smaller branching fraction, at least 
less than the OZI allowed processes, e.g. J/ip — ► j<p<f> or -fuiu. It seems that the X(1812) has a 
large contribution to the J/ip — ► joj(f> decay. Moreover, if the scalar of X (1812) is assigned as the 
qq ordinary meson, it weakly couples to the decay mode of uj<p and has a small phase space to decay 
into u>4> near the threshold. However, the observed branching fraction is at the order of 10~ 4 . So 
it is difficult to fit the X(1812) into the spectrum of ordinary qq mesons. Many efforts have been 
made to interpret the X(1812) as an exotic state [21 [3J H], such as glueball|5l [3], hybrid [71 [5] and 
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four-quark state [9]. The contributions from final state interactions of /o(1710) decays into PP 
and VV states are already investigated [TO] . 

Search for ordinary decay modes of the X(1812) is essential for understanding its structure. Like 
the other scalar mesons, e.g. /o(1500), one expects that the Jf(1812) may have a large fraction 
to decay into PP and VV states. However, no more experimental information on AT(1812) are 
available, a thorough investigation on all possible production and decay mechanisms for AT(1812) 
is thus necessary for understanding data, such as the molecule scenario should be inspected. The 
wide resonance K* and its pair (K*K*) near the A(1812) threshold make it a candidate for such 
consideration. PDG [H] quotes M = 891.66 ± 0.26 MeV and T = 50.8 ± 0.9 MeV as the average 
value of the mass and total width for the K* ± 1 and almost the same values for the K*°. If the 
width of a resonance is regarded as the full width at half maximum of the mass distribution in 
a nonrelativistic Breit-Wigner form, the K*K* pair still has a large probability to lay within the 
V(1812) mass region. Moreover, the production of K*K* is copious in the J/tp — > r yK*K* with 
the branching fraction Br(J/ip — ► ^yK*K*) = (4.0 ± 1.3) x 10 -3 [16]. Some interactions between 
the K* and K* can not be ruled out, thus give rise to the reaction K*K* — > ui(f) possibly. 

The theoretical studies on the scalar meson structure are always the controversial subjects in 
particle physics, and the sparse experimental information makes it more difficult to draw a decisive 
conclusion about their nature. Even for the well established scalar / (980) and a (980), there are 
still more puzzles about their masses, decay rates and so on. Many theoretical efforts have been 
made in literatures to study the scalar as the exotic meson, such as the four-quark state [T7] , glueball 
or hybrid [TS] and molecular state [15]. In the quark potential model, the molecule of K*K* is 
predicted to be /o(1710) [20] with a bound energy about 70 MeV. However, the quark model based 
on the pairwise effective interactions predicted that there exists a weak bound molecule near the 
threshold, such as the DD molecule [21] with the bound energy about 3 MeV. Near the K*K* 
threshold, the possibility that some interactions between K*K* pair might give rise to a weakly 
bound molecular state can not be ruled out. 

In brief, the production and decay mechanism of the A"(1812) observed in the J/tp — > jui(p 
deserve to do a thorough investigation, especially by assigning the A"(1812) as a molecular state, 
from which we expect to learn more about X (1812). The (K *K*) pair could be a candidate of such 
consideration. As follows, we will formulate the AT(1812) decay and production rates by assuming 
the X(1812) as the K*K* molecule, from which we hope to see whether the experimental data 
could be understood in this scenario, and which channel could be possibly used to look for the 
V(1812) in experiment and to find its trace again. Therefore, all possible descriptions about the 
V(1812) structure are subject to the experimental test in the future, and the comparison between 
the data and the theoretical prediction may shed some light on the X(1812) nature. 

2 Decay Mechanism 

The experimental information on the X(1812) shows that the J PC = ++ is favored. So, in the 
scenario of (K*K*) molecular state, the K*K* is chosen as the orbital ground state L = and 
the spin singlet S = due to the restriction of the space-orbit parity [P — (— 1) L ] and the charge 
parity [C — {— l) L+s ], the wave function is expressed as [213] : 

*(K*K*) = C, lmi,lm2 e ( TO l) e ( TO 2)^*(Pl)^-(P2)*(X*i?«)(Pl -P2), (1) 

mi ,7712 

where e(rrii) is the spin wave function of K* or K* , and the (pn(pi) is the K* or K* wave function 
in momentum space. 

At hadronic level, the K*K* can be scattered into the u>(f> final sate by exchanging a meson 
as assumed in [10] . At quark level, such a meson- meson scatterring process has already been 
investigated by many authors in quark model [ll]-[15j. It is found that the process can be described 
by one gluon exchange (OGE) between qq or qq pair, by which one found that it gives an excellent 
description both for the light and heavy meson-meson scattering process. By analogy with these 
models, we describe the reaction K*K* — > MM to begin with the quark-quark scattering, then 
leads to the quark rearrangement between the two color singlet clusters and subsequent formation 
of the final two mesons as shown in Fig. [IJa-d). In addition, the qq annihilation may take place 
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Fig. 1: The four meson-meson scattering diagrams by exchanging f\fc quarks, 
and similarly by exchanging f\f\ quarks. 
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Fig. 2: The schematic diagram for decays of the K*K* molecule via a qq annihilation. 



between K* and K* clusters and leads to the (K*K*) molecule to decay into the different final 
state as shown in Fig. [2] 

As the K*K* in the molecule are weakly bound, we may relate this decay rate to the cross 
section a(K*K* — ► M\M 2 ) near threshold, by analogy with the calculation of the decay rate of 
So positronium to 77. The decay rate for a (i4T*^*)-molecule is |22) : 

T((K*K*) -> MiM 2 ) = a(K*K* -> M X M 2 ) ■ v rel ■ \^ (K , Rn {0)\ 2 , (2) 

where v re i is the relative velocity of the K* and K* , ^ik'K*)(^) 1S tne (K* K*)-mo\ecu\e wave- 
function at the origin r = 0. The cross section can be expressed by the invariant amplitudes for 
the scattering of K*K* to the final meson pair, i.e. 

M(K*K*) -> M,M 2 ) _ 1 1 IA „^™ ,,,,2 



dt 647TS 2 \P A cm\ 2 



\M((K*K*) — > MiM 2 )\ , (3) 



where the Mandlestam variable t — M\ — 2EaE c + 2P A • P c + Mq. The invariant amplitude M. 
is generally expressed as: 

M = (M 1 M 2 \H I \(K*K*)) = F color ^ avor ^ spin+space (4) 

For the quark-exchange mechanism and the quark-annihilation mechanism, on has F co j or = 4/9. 

The flavor factors are calculated based on the flavor functions of the final meson pairs with the 
interaction Hf, which can be written as, for example, 

„ _ / &hhh*U S fih S hh for Fig. [1] (a), 
1 \ «; ; " ;."/.«./, for Fig. 1. w 

where / denotes quark flavor, and aj are the quark / creation and annihilation operator, 
respectively. 

The contributions from the spin and space parts are denoted by: 

-^spin+space = ( i PM 1 fM 2 \H ss \tf K ,Lp R ,), (6) 

where </?m and tpx* denote the wave functions of the decayed mesons and K* molecular state in spin 
and coordinate space, respectively. They are directly evaluated by Feynman diagrams as shown in 



Fig. [T]and Fig. [2] It reads 
1 
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For the quark /1/2 exchange mode (Fig. [1]), the operator Ti reads: 
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where u(si,Pi) and v(si,pi) are the quark and antiquark Dirac spinor with the normalization 
condition u' (s,p)u(s,p) = —v' i (s,p)v(s,p) = 2p°. For the quark /1/2 exchange mode, it has the 
same diagram and the operator Ti ex takes the similar form. 

For the quark /1/2 and /1/2 annihilation mode (Fig. [J), the interaction Ti are respectively given 

by: 

;v(s 7 ,p 7 )-f ll u(s 6 ,p 6 )u(s 3 ,p 3 )-f> Jl v(s2,P2)Ss u sJs 4 ,s 5 5 3 (p 1 - p 8 )<5 3 (p 4 ~ Ps), (12) 
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3 Production Mechanism 

The production of the molecule (K*K*) is assumed via the J ftp radiative decay into a real photon 
(7) plus two virtual gluons (gg), followed by the formation of the meson pair of K*K* , and then 
possible formation of the K*K* molecule through the final state interactions occurs. To note 
that the decay J/tp — > jK*K* has a large branching fraction (Br=(4.0 ± 1.3) x 10 -3 ), then the 
interactions between the K*K* pair may take place and may give birth to the (K*K*) molecule. 




Fig. 3: The schematic diagram for the production of the 
(K*K*) molecule in J/ip radiative decays. 



To evaluate the production rate of (K* K*)-molecule, we consider the process as shown in Fig. [3] 
at the level of the leading order of perturbative QCD, and the bound states are phenomenologically 
described by their wavefunctions. The whole amplitude can be decomposed into two parts as done 
in 1231 [Ml , which is written as 



K l K 2 



(14) 



where k% and &2 are the momentum of the two virtual off-shell gluons. The matrix elements of 
(J/tp\igg) an( i (gg\(K*K*)) describes the subprocess J/ip — > -fgg and gg — > (K*K*), respectively. 
The sum is over the polarization vectors of the two gluons. 



For the subprocess J/ip(E) — ► 7(e, k)g(ei, ki)g(e2, k 2 ), the amplitude can be expressed by [23] 
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with 



o~u = { e i • e 2l — ^1 ' ^ e ' • ki — k 2 ■ fee* • k\E ■ k 2 — k\ ■ kk 2 ■ kE ■ e*] 
+ E ■ e*[k\ ■ ke\ ■ k 2 e 2 ■ k + k 2 ■ ke 2 ■ kit\ ■ k — k% ■ k 2 e\ ■ ke 2 ■ k]} 
+ {ex,ki <-> e,k} + {e 2 , k 2 e, k}, 



(16) 



where E, e\ and e 2 are polarization vectors for J/ip, the real photon and the two virtual gluons, 
respectively, and their momentum arc denoted by k, k\ and k 2 , respectively. R^(0) is the J/ip 
radial wave function at origin in coordinate space. 

For the subprocess g(e\, ki)g(e 2 , k 2 ) — > K* (ek» , P 2 )K* (e^, , P 2 ), the amplitude can be directly 
calculated from the leading-order diagram as shown in Fig. [3] by using the standard Feymman 
rules. The details of the calculation are given in appendix A. The amplitude is calculated to be: 
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where R R * (0) is the K* wave function at origin, and p\ and q\ are respectively the momentum of 
u— quark and s— quark bound in K* , and the antiparticle counterparts denoted by p 2 and q 2 with 
Pi = m u Pi/(m s + m u ) and q t = m s Pi/ (m s + m u ). 

If we assume the K*K* pair to be a molecular state, which is phenomenologically described by 
a wavefunction <5(Pi — P2) given in Eq. (fTJ, then the amplitude for gg — ► (K*K*) reads: 
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(18) 



After inserting Eq. (fT5"|) and (jXTJ) into Eq. (JT^J) , and making substitution k\ — p\ + q 2 and 

k 2 = p 2 + qi, one obtains the amplitude for J/ip — > jK*K*. Using the standard formula for 
three-body decays [IB], the decay width for J/ip jK*K* is written as: 
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where 03 is a standard 3-body phase space factor. 

Similarly, the decay width for J/ip — > j(K*K*) is written as: 
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where P 7 is the photon momentum. 



4 Numerical Results 

To evaluate the (K*K*) molecule decay width and the production rate, one should have a reliable 
method to describe the wavefunctions of the low lying mesons. However, no rigorous theory from 
the QCD first principle is available to describe the light bound states, so one has to use a phe- 
nomenological model to include the nonperturbative properties. Phenomenologically, we construct 
meson wave functions in the constituent quark model; they are decomposed into three parts, i.e. 
the flavor, spin and the space wavefunctions. For example, the wave function of K* + is constructed 
as: 

\K* + ) = \us)e(l,m)ip(p), (21) 



Tab. 1: The flavor factor for the K*K* molecule decays into 
a meson pair via the quark exchange and annihilation mech- 
anism, 
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Tab. 2: The harmonic parameters calculated with the mass 
M and the decay constant fu [26] by the relation (5 = 

V^(i) 1/3 /M /3 
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where e(l, m) denotes the symmetry spin wave function of two constituent quarks. For pseudoscalar 
mesons, it is chosen as asymmetric spin wave functions. 

For rj and rj' mesons, we choose the mixing scheme to expand their flavor functions in terms of 
the singlet and octet quark flavor basis as: 

|r/) = cos0|t7,j) - sin0|Tk), 

= sin0|ry g ) + cos0|?7 s ), (22) 

where is the mixing angle, \rj q ) = \uu+dd) j v2y 9 (p) and \r) s ) — \ss)<p a (p), where y(p) is chosen as 
the ground harmonic oscillator basis associated with the decay constant f q or f s . In [5S] the mixing 
angle and the decay constants are extracted from the experimental data as: = 39.3° ± 1.0°, f q = 
(1.07±0.02)/ OT and f s = (1.34 ± 0.06)/*. The flavor factors for X(1812) rtf, r]r], together with 
the other final states, are given in Table [T] 

The space wave function ip(p) is solely dependent on the harmonic oscillator parameter (3. From 
some studies on the the spectroscopy and decay rates of the low lying scalar and vector mesons, 
their decay constant is determined; they are related to the meson wavefunction at the origin by 
fp.v = \J Mp V I0p,v(O)|- On the other hand, the wave function at the origin is related to the 

parameter by F ,y(O) = (/3 2 /vr) 3 / 4 . Table OH summarizes the decay constants determined by 
experimental data 26] and the harmonic oscillator parameters determined by the relation f3p y = 

For X(1812) — > M1M2 decays, the decay width includes the contributions from the quark 



Tab. 3: The production rate of F^t K *x^/T in terms of 

the rms of the (K*K*) molecule. 
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exchange and quark annihilation mechanisms, respectively denoted by T|^ M and T^ M2 . Due to 
the non-zero flavor factors, the quark exchange mechanism allows only three decays, uxj), r)T) and 
7777' channels as shown in Table [TJ The calculation of the decay width is straightforward though 
tedious according to Eq. The numerical calculation yields the ratios: 



pex . pex . pex . pex . pex . pex . pex . pex _ 

1 r/n' ■ 1 rin ■ 1 Lj<f> ■ 1 K+K- ■ 1 p+p- ' 1 lolo ■ 1 K' + K'- • 1 tt+tt- — 

9.3 : 4.3 : 1.0 : : : : : 0. (23) 

This indicates that the quark exchange mechanism enhances the decays X(1812) -+ r\r\ and 7777'. On 
the other hand, the quark annihilation mechanism makes no contribution to the decay X(1812) — > 
uxj). To determine the contributions from this mechanism, we estimate the reduced ratios r an 
T a ^ M2 /r^. The numerical calculation yields: 



MiM 2 



pan . pan . pan . pan . pan . pan . pan . pan _ 

1 m , ■ i- nv ■ 1 U( p ■ 1 K+K- ■ 1 p+p- uju ' K' + K'- ■ 1 -K+-K- — 

0.39 : 0.04 : : 0.44 : 0.41 : 0.005 : 0.001 : 1.9 x 10~ 5 . (24) 

This indicates that the quark annihilation mechanism makes small contributions to X(1812) — > 
M1M2, especially for the 7r + 7r~ mode, which is highly suppressed in dynamics. If the both decay 
mechanisms and the interference effects between them are taken into consideration, the total decay 
widths for each mode are calculated to be: 



r?jtj' : : : r K + K - ■. r p + p - : r ww : Tk*+k*~ 



13.6 : 5.14 : 1 : 0.44 : 0.41 : 0.005 : 0.001 : 1.9 x 10~° (25) 

For the production of the (K*K*) molecule in the J/tp — * j(K*K*), the production rate depends 
on the overlap of the K* and K* wave functions in momentum space. They are phcnomcnologically 
described by ground harmonic oscillator bases with a parameter /3, which is related to the root mean 
square (rms) radius (r 2 ) by the relation of (3 = . The model study using a nonrelativistic 

Hamiltonian with pairwise effective interactions shows that the rms becomes large when the binding 
energy of two- meson molecule decreases [21] , With the binding energy ranging from 100 MeV to a 
few Mev, the rms runs from 0.3fm to 2 fm. Near the threshold of two bound mesons, the rms of 
the molecule is about 1-2 fm. In our calculation, we naively set the (3 parameter within the range 
0.5 ~ 2.0 fm. From Eq. (|19H20p the numerical calculation yields the ratios of the decay width for 
the J/ip — > ~fK*K* and j(K*K*), which are given in Table The results show that the strong 
binding of the (K*K *) molecule favors its production from J/tp decays. 



5 Discussion and Summary 

Based on the scenario of X(1812) to be a (K*K*) molecule, the decay width and the production 
rate are evaluated. It is found that the quark exchange mechanism plays an important role in the 
understanding of the decay X(1812) — * uj<j>, and we find that this decay mechanism enhances the 
decay widths for 7777 and 7777' modes. On the other hand, the quark annihilation mechanism makes 
a small contributions to X(1812) — > PP and VV decays, and the decay X(1812) — * tttt is highly 
suppressed in dynamics. These results suggest that if X(1812) is a (K*K*) molecule, it should 
be found in the decays J/ip — > jX, X — > 7777' and 7777 except for the observed mode X — > uxj). 
Unfortunately, due to the low statistics of the present data sample, no information on J/ip — > 77777' 
and 77777 decays are available in PDG table. 



Figure [4] shows experimental information on the invariant mass distribution of two mesons in 
J/ip — ► 'ycocp, JTjrj, "fp + p~, jK + K~, "ftoLu, ~/K* + K*~ and 77r + 7r~ . In X{ 18 12) mass region as seen 
in the J/ip — > 7w</>, no signal evidence is observed in the other decays. As for the decay J/ip — ► 77777, 
it is firstly reported by the collaboration of Crystal Ball detector twenty years ago with very low 
statistics, but after that no confirmation was made by other collaborations. Our predictions deserve 
to be tested in this decay, together with the J/ip — > 777?/. With a large data sample accumulation 
and the improvement in the detector performance, especially for the photon identification for BESIII 
and CLEOc detector, search for X(1812) signals in these two channels are possibly achieved. In 
the J/ip — > 7/9 + p~, the partial wave analysis (PWA) shows that the contribution from scalar 
near the mass of 1690 MeV dominates the m pp mass distribution. In the J/tp — > ^tt + tt~, the 
PWA shows that significant resonance with mass 1270 Mev favors the J PC = 2 ++ , and small 
scalar components with mass equal to 1466 MeV and 1765 MeV are also reported. Recently, the 
BESII collaboration reported the PWA performed on the J/ip — ► ^uslo. The results show that the 
contributions in the mass distribution of m^u below 2.0 GeV are predominantly from pseudoscalars, 
and only with small contributions from / (1710), / 2 (1640) and / 2 (1910). The PWA of the decay 
J/ip — > r yK + K~ is also reported by the BESII collaboration, the dominant contributions are from 
the scalar / (1710) and the tensor /£(1525). The PWA on the decay J/ip -> jK*K* is reported 
by BESI collaboration with 7.8 million J/ip data, the results showed that a broad h resonance 
with mass M = 1800 MeV is observed, no significant ++ , 1 ++ or 4 ++ signal is found. In short, no 
evidence for X(1812) is currently observed in the hadron mass distributions of the radiative decays: 
J/ip — > 7/9 + p _ , r yK + K~ , jojuj, 'fK* + K*~ and 77r + 7r~ , from which it seems to be consistent with 
the our calculation based on the (K*K*) molecular picture. The further test of this scenario is 
expected to search for the X(1812) in J/ip — > 77777' and 77777 decays in the future. 

If we naively assume that the [K* + K*~) molecule dominantly decays into 7777', r)r),u>(p, KK 
and pp, then one obtains Br(X — * ucp) ~ 4.60% from our calculations. From the measurement 
of Br(J/tp -> jX)Br(X -> ucp) = (2.61 ± 0.71) x 10~ 4 , one has the production rate Br(J/ip -> 
7V) = (5.67 ± 0.07) x 10~ 3 . Combined with the PDG value, one has T(J/ip -> jX)/T(J/ip -> 
~fK*+K*-) = 2.83±0.92. Our calculation yields T(J/ip -> jX)/T(J/ip -> jK*+K*-) = 2.13±J;||, 
where the central value corresponds to the \J (r 2 ) =1.0 fm, and the uncertainty corresponds to the 
range 0.5 fm < J (r 2 ) < 2.0 fm. 

It is also important to look for decay rates of X(1812) — > PP and VV in experiment for 
understanding its nature. Theoretically, these decays have been evaluated numerically by many 
authors using different models to describe the X(1812) structure. For example, in four-quark 
picture [S], it is found that V(1812) — > uolo, K*K* are the two dominant decay channels and 
X(1812) — > KK, r\r\, ryrj are suppressed. While in glueball picture [5], the dominant decays are 
X(1812) — * pp and u>u>, and decays into tocp and K*K* are highly suppressed. In the quarkonia- 
glueball-hybrid mixing scheme [5], it turns out that the branching fractions for X(1812) — > KK 
and 7777 are about 30%, which adds up to about 70% of the total V(1812) decay width. However, 
in our (K*K*) molecule picture, the dominate decays are X(1812) — > 7777' and 7777. 

In summary, we evaluate the decay widths of X(1812) —> 7777', 7777, uj<p, K + K~, p + p~, ujuj, K* + K* 
and 7r + 7r _ based on the assumption of the JT(1812) as a candidate of (K*K*) molecule. It turns 
out that the quark exchange mechanism plays an important role in the understanding of the large 
decay with for the X(1812) —> ucp. We also find that the decay widths for V(1812) — > 7777' and 
7777 decays are enhanced by the quark exchange mechanism. These channels are suggested to be 
the tools to test the molecular scenario in experiment. The branching fraction of Br(X — > ujcp) is 
evaluated to be about 4.60%. Searches for additional evidence about the V(1812) in J/ip radiative 
decays are reviewed. In the molecular scenario, the X(1812) production rate is also evaluated to be 
T{J/ip -> jX)/T(J/ip -> 1 K*+K*-) = 2.13+^, which is close to the measured value 2.83 ± 0.92. 
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Fig. 4: The experimental measurement of mass distributions m/ ll / l2 for 
radiative decays of J/ip — > 7/11^2, (a) <^A/> - ► 7k><ft[27]. (b) J/V' — > 7W 
[28], (c) J/^ -> 7P + p" [29], (d) J/?/> -» jK + K~ [30], (e) J/^ -> 
7ww^l, (f) J/V> -> 7^*+^*" [32], (g) J/^ -> 77T+7T- [33]. 



A The amplitude for gg — ► K*K* 

We start with the process g(e 1 )g(e 2 ) -> g(pi, s x )q(p 2 , s 2 )q(q 2 , S2) -» # *(Pi)^*(P 2 ) as 

shown in Fig. [31 where (i = 1,2) denotes the polarization vector of the gluon; the momentum 
and the spin of the quark (antiquark) q (q) are denoted by Pi (qi) and (sj), respectively; the 
outgoing momentum for K* (K*) is denoted by P 1 (P 2 ). The amplitude with Lorentz indexes Hi 
and /i2 is defined as: 

1 d 3 pi d 3 q 1 d 3 p 2 d 3 q 2 



3^ K * K * J J2p\2q«2pppi (2tt)3 (2tt) 3 (2nf (2ir) 3 

x ^ u(s 1 ,p 1 )-y f * 1 v(s 2 ,q 2 )u(s 2 ,p 2 )-f> 12 v(s 1 ,q 1 ) 

all spin indexes 

X ^ Si; \ s i\SiS lz )(^s 2 ; ^s 2 \S 2 S 2z )(L 1 M 1 ; SxSu\JiJu)(L 2 M 2 ; S 2 S 2z \J 2 J 2z ) 

X */<*L 1 5 1 (Pl-qi)*i?.L 2 S 2 (P2-q2) 

x (2 7 r) 3 (5(P 1 -p 1 -q 1 )(2^) 3 ( 5(P 2 -p 2 -q 2 ), (26) 

where ^k»l 1 s 1 k*l 2 s 2 ) denotes the wave function for K* (K*) with the orbit and spin angular- 
momentum quantum number L\ (L 2 ) and Si (ft), respectively. For the vector meson K*, it is 
assigned as the 3 Si state in qq quark model. So we have Li = and Si = Jj z = 1 (i = 1, 2). 

The above amplitude can be simplified with the help of the spin projection operator Pslsl z -s 2 S 2z > 
which is defined as: 

P s'sLs 2 s 2z = E "(suPih^vfr, d2)u{s 2 ,p 2 )r 2 v{s uqi ) 

s 1 s 1 s 2 s 2 

x (j s i; ^s 1 \S 1 S lz )(^s 2 ; ^s 2 \S 2 S 2z ) 



= Tr 



^2 w(si,9i)u(si,pi)7 Atl (isi; ^-s 1 \S 1 S l2 



x ^2 v ( s 2, q2)u(s 2 ,p 2 )j^ 2 (-s 2 ; -s 2 \S 2 S 2z 

s 2 s 2 

= Tr[P 1Sl _,P 1S2 J, (27) 

where Pis lz and Pis 2z are the spin projection operators for the spin-1 particles. As given in [34], 
they are explicitly expressed with the spin polarization vectors 6k * as: 

p iSz(Pi><li) = ^ v {si,qi)u{s l ,p i )Y{]^Si^s l \lS z ) 

Si Si 

1 

(-ii + m s )e K *(S z )(j!> + m u ). (28) 



y/(qf + m s ){pf + m 

To simplify the integral in Eq. ([2"6"[) . one method to make an approximation is the nonrelativistic 
assumption for valence quarks, ie. pi — qi = 2q — * 0. Then the integral of the amplitude can be 
evaluated at the origin with the substitution of the momentum n, = m ™ Pi and o, = ma , Pi , ie. 



f d 3 p! d 3 qi 3 Rk*o(0) , on , 

7 ( 27 r)3 { 2n)^ K ' 0s{Pl qi)(2?r) * (Pl Pl qi) * (29) 

where i?^*o(0) is the radius wavefunction for the K* at origin. With the help of equations ([2"7]) to 
the amplitude for gg — > K*K* can be simplified as: 

2! VpSp§9?9§ + ™s)(<7 2 ° + m fl )(p? + m u )(p° + m„) 4^ 

x Tr[(-^i +m s )^K*{i>\ +m u )£ 2 (-q ; 2 +m s )e R ,(p' 2 +m u )/ti 
+ (ei^e 2 )], (30) 

where the factor of 1/2! comes from the contribution of the cross term for the two gluons of identical 
particles. 
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